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Abstract 

We show the existence of an energetic solution to a quasistatic 
evolutionary model of shape memory alloys. Elastic behavior of each 
material phase/variant is described by polyconvex energy density. Ad¬ 
ditionally, to every phase boundary, there is an interface-polyconvex 
energy assigned, introduced by M. Silhavy in [44]. The model con¬ 
siders internal variables describing the evolving spatial arrangement of 
the material phases and a deformation mapping with its first-order gra¬ 
dients. It allows for injectivity and orientation-preservation of defor¬ 
mations. Moreover, the resulting material microstructures have finite 
length scales. 

Keywords: poly convexity; shape memory materials; rate-independent 
problems 


1 Introduction 

In elasticity theory, it is assumed that experimentally observed patterns are 
minimizers or stable states of some energy. Shape memory alloys in particu¬ 
lar have a preferred high-temperature lattice structure called austenite and 
a preferred low-temperature lattice structure called martensite. Such shape 
memory alloys, as e.g. Ni-Ti, Cu-Al-Ni or In-Th, have various technological 
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applications, for an overview see e.g. [MJ. The austenitic phase has only one 
phase/variant but the martensitic phase exists in many symmetry related 
phases/variants; the mixing of these different phases can lead to the forma¬ 
tion of complex microstructure. In the continuum theory, the total energy 
of the system is described in terms of a bulk energy which describes elastic 
stresses and an interfacial energy, concentrated on the interfaces between 
the different phases. We establish existence of quasistationary solutions for 
a model, where it is assumed that the bulk part of the energy is polyconvex 
while the interfacial part of the energy satisfies a corresponding condition of 
interfacial polyconvexity introduced by Silhavy [441 145] . The model describes 
the evolving spatial arrangement of the material phases and the deforma¬ 
tion of the sample. It allows for injectivity and orientation-preservation of 
deformations. Moreover, the resulting material microstructures have finite 
length scales. 

To investigate the existence of a global minimizer of the energy for static 
variational problems from elasticity, different notions of convexity have been 
considered. For problems with a single material phase, a well justified no¬ 
tion of convexity which is sufficient to ensure the existence of a minimizer 
is the notion of polyconvexity due to Ball El El- It is also fairly easy to 
construct examples of polyconvex functions which makes it attractive for 
continuum mechanics of solids. On the other hand, in shape memory alloys, 
many different phases might coexist. If interfacial energy is not taken into 
account, then global minimizers of the energy in general do not exist. A 
way out is to use relaxation methods, searching for the so-called quasicon- 
vex envelope of the specific stored energy mm or using Young measures 
[25 , 26l f3lj 38] . Let us point out some partial results which have been 
obtained in this direction: We refer to jTj for a weak* lower semicontinuity 
results for sequences of bi-Lipschitz orientation-preserving maps in the plane 
and to [6] for an analogous result along sequences of quasiconformal maps. 
Then [32] found relaxation including orientation preservation for p < d, 
where d is a spatial dimension. Finally, [HI derived a relaxation result for 
orientation preserving deformations with an extra assumption on the result¬ 
ing functional, namely that the quasiconvex envelope is polyconvex. There 
also exist various phenomenological models of shape memory alloys which 
are convenient for numerical computations; see e.g. [8j. 

On the other hand, models have been considered where interfacial energy 
is taken into account. Such models have been e.g. used to estimate the 
scaling of the minimal energy and to derive typical length scales of patterns. 
The minimal scaling of the energy of an austenite-martensite interface has 
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been studied by Kohn & Muller and Conti in [30j [31, 16] for a 2-d model 
problem, the three-dimensional case and more realistic models have been 
investigated e.g. in Ha m mm m m, for similar analysis on related 
models see e.g. mmm- In these models, either a .BC-penalization of the 
interfacial has been used or a penalization of some L p -norm for the Hessian 
of the deformation function. In general, the specific form of the energy is, 
however, not clear from physical considerations. In the literature, necessary 
and sufficient conditions for the specific form of the interfacial energy have 
been investigated recently which allow for the existence of minimizers [2lL 
02]. Recently, Silhavy has introduced a notion of interface polyconvexity and 
has proved that this notion is sufficient to ensure existence of minimizers for 
the corresponding static problem @31145] , In this note, we extend this static 
model to a rate-independent evolutionary model and prove existence of an 
energetic solution. 

In shape memory allows, the stored energy density W : M 3x3 — > M is min¬ 
imized on wells SO(3)T), i = 0,..., M , defined by M positive definite and 
symmetric matrices Fq,, Fm , each corresponding austenite and M vari¬ 
ants of martensite, respectively. By the choice of reference configuration, we 
may furthermore assume Fq := Id (the identity), i.e. the stress-free strain of 
austenite is described just by the special orthogonal group SO(3). In non¬ 
linear elasticity, the energy density W is usually formulated as a function 
of the right Cauchy-Green strain tensor F T F. Note that this tensor maps 
the whole group 0(3) of orthogonal matrices with determinant ±1 onto the 
same point. Thus, for example, F H > \F T F — Id | is minimized on two en¬ 
ergy wells, i.e., on SO(3) and also on 0(3)\S0(3). However, the latter set is 
not acceptable in elasticity since corresponding deformations do not preserve 
the orientation. Additionally, notice that, for example, considering arbitrary 
Q £ 0(3) \ SO(3) and an arbitrary R £ SO(3) such that Q and R are rota¬ 
tions around the same axis of the Cartesian system then rank(Q — R) = 1, 
i.e. Q and R are rank-one connected and determinant changes its sign on 
the line segment [Q',R\. Convex combinations of rank-one connected ma¬ 
trices play a key role in relaxation approaches of the variational calculus 
@] [5] DEI [33]. This shows that is it important but also not straightforward 
to ensure that solutions of static or evolutionary problems are physically 
sound, in the sense that they preserve orientation. This is, in particular, 
unclear on models based on quasiconvexification as described above, since 
usually there is no closed formula of the envelope at disposal and since 
physically justified conditions on deformations as orientation-preservation 
and injectivity are not included in these models. On the other hand, our 
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solutions are constructed in a way such that the obtained time dependent 
deformations are orientation preserving and injective and no additional reg¬ 
ularization of variables is needed if passing from a static to an evolutionary 
model. Injectivity and orientation preservation is not very often considered 
in the theory of rate-independent processes. We refer to [37] for treatment 
of a model in nonlinear elastoplasticity. 

Structure of the paper: After introducing used notations, in Section 2, we 
first describe our model, the stored elastic energy, loading, and dissipation. 
In Section 3, we state and proof our main result, the existence of an energetic 
solution. As it is nowadays a standard procedure; cf. e.g., [22] we only sketch 
the main steps and pay more attention to injectivity of deformations which 
is not frequently treated in the frameworks of rate-independent evolutions. 
We refer, however, to [37] for numerical approaches to finite elastoplasticity 
including injectivity. 

Notation: The spaces 4F 1,P , 1 < p < oo, denote the standard Sobolev 
space of L p -functions with weak derivative in L p . Further, BV stands for 
the space of integrable maps with bounded variations, see e.g. [1] [20] for 
references. For a (measurable) set E Cl 3 , we denote its three-dimensional 
Lebesgue measure by £ 3 (E) and its two-dimensional Hausdorff measure by 
7i 2 (E). The space of vector valued Radon measures on hi with values in Y 
is denoted by M(Cl, Y). 

Let fi C O C K 3 be Lebesgue measurable sets and let B(x,r) := {a £ 
R 3 : \x — a| < r} . For x £ 0 we denote the the density of Cl at x by 
9(Cl,x) := lim. r _>.o C 5 (Cl H B(x, r))/C^(B(x, r)) whenever this limit exists. A 
point x £ is called point of density of fl if 9(Cl,x) = 1. If 9(Cl,x) = 0 
for some x £ 0, then x is called point of rarefaction of hi. The measure- 
theoretic boundary d*Q of 0 is the set of all points x £ hi such that either 
9(Cl,x) does not exist or 9(Cl,x) fL {0,1}. We call Cl a set of finite perimeter 
if E 2 (d*Cl) < Too. Let n £ R 3 be a unit vector and let H(x,n ) := {x £ 
0 : (x — x) ■ n < 0}. We say that n is the (outer) measure-theoretic normal 
to Cl at x if 9(Cl I~1 H(x, —n),x) = 0 and 0((fi \ Cl) n H(x , n),x) = 0. The 
measure-theoretic normal exists for B 2 almost every point in <9*0, see e.g. 
[201116]. 

For two matrices A = ( aij),B = ( bij ) £ R 3x3 , we define A : B = atjbij 
with Einstein’s sum convention. By Ixnwe denote the tensor defined by 
(A x n)b = A(n x b ), i.e. (A x n)kj = where euj is the Levi-Civita 

symbol. One can easily check that the cofactor matrix of A £ R 3x3 in 
terms of the Levi-Civita can be expressed as cof A = ^ (eikN-jpqO-kp a (:q)ij■ In 
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particular, we get d (lkt (cof A)^ = \d a . kf {e ik e^jp q a k pa eq )ij = eik q eji p a qp . We 
refer e.g. to [231 for a definition of the surface gradients V 5 . If n G M 3 is 
an outer unit normal to the surface S , then V 5 := V(Id—n®n), where we 
recall that Id denotes the unit matrix in M 3x3 . 


2 Model description 

2.1 Elastic energy 

Admissible States: We assume that the specimen in its reference configu¬ 
ration is represented by a bounded Lipschitz domain QcM 3 . We consider 
a shape memory alloy which allows for M different variants of martensite. 
The region occupied by the i-th variant of martensite is described by the set 
Oj C 0 for 1 < i < M, while the region occupied by austenite is given by 
Po C II. We assume that the sets f2j are open and have finite perimeter. Fur¬ 
thermore, the sets are pairwise disjoint for 0 < i < M and N := P\(J. Pj 
is a set of zero Lebesgue measure. The case 0* = 0 for some 0 < i < M is 
not excluded. The partition of Q into {H-di'io can be then identified with 
a mapping 2 : fl —> R M+1 such that Zi(x) = 1 if x G and Zi{x) = 0 else. 
We call z the partition map corresponding to {Qi}f £ 0 - Clearly, with the sets 
Pi chosen as before, we have 0 z i( x ) = 1 f° r almost every i£fl and the 

function z is of bounded variation. We hence consider z £ Z, where 


M 

Z := jz G BV(P, {0,1} M+1 ) : ZiZj = 0 for i / j, ^ z t = 1 a.e. in j. 

i=0 

In order to describe the state of the elastic material, we also need to intro¬ 
duce the deformation function y G IT 1 ,P (P,M 3 ), p > 3, which describes the 
deformation of the elastic body with respect to the reference configuration 
P. We hence consider deformations y £ y, where 


y = jy € W 1,P (J1,R 3 ) : det Vy > 0 a.e. , f det'Vy(x)dx< 


where we will always use the assumption p > 3. The integral inequality to¬ 
gether with the orientation-preservation is the so-called Ciarlet-Necas condi¬ 
tion which ensures invertibility of y almost everywhere in P PUCES]. In the 
following, we will assign to each state of the material (y,z) G y x Z an elas¬ 
tic energy £. In our model, the energy consists of a bulk part Ey,. penalizing 
deformation within the single phases, an interfacial energy E- mt , measuring 
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deformation of the interfaces between the phases and a contribution L(t,-) 
which measures work of external loads, i.e. 


£{t,V,z) ■= E h (y, z) + E int (y, z) - L{t,y). (1) 

Here t denotes time to indicate that we will deal with time-dependent prob¬ 
lems. We will specify these three parts of the energy in the following. 

Bulk energy: The total bulk energy of the specimen has the form 

Eb(y,z):= [ W(z(x),Vy(x)) dx, ( 2 ) 

Jn 

where we assume that the specific energy W : R A/+1 x R 3x3 —y R U {+ 00 } 
of the specimen can be written as 

M 

W{z,F):=Y,ZiWi(F)=:z-W{F), (3) 

i=0 


where W{, 0 < i < M, is the specific energy related to the z-th phase of 
the material and W := (Wo,...,, Wm)- We will work in the framework of 
hyperelasticity, where the first Piola-Kirchhoff stress tensors of austenite and 
martensite have polyconvex potentials denoted by Wo (austenite) and Wi, 
i = 1 ,,M for each variant of martensite, see e.g. [Mj and the references 
therein. For 0 < i < M, we therefore assume 


Wi(F) : = 


hi ( F , cof F, det F) 
+00 otherwise 


if det F > 0, 


( 4 ) 


for some convex functions hi : R 19 —> R. We use the following additional 
standard assumptions on the specific bulk energies W t . For 0 < i < M and 
F e r 3x3 , we assume that for some C > 0 and p > 3 

Wi(F) >C(-1 + \F\ P ) VFe R 3x3 , (5) 

Wi{RF) = Wi{F) Vi? 6 SO(3), F e R 3x3 , (6) 

lim Wi(F) = + 00 . (7) 

det F—¥ 0+ 


Interfacial energy: We consider the interfacial energy in the form introduced 
by Silhavy in HU S5J: We hence assume that the specific interfacial energy 
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fij between the two different phases i,j £ {0,..., M} can be written in the 
form 

n) = gi(F, n ) + gj(F, n ), (8) 

where F £ R 3x3 and n £ R ' 3 is a unit vector such that Fn = 0, We assume 


gi(F,n) := i(n,F x n,cofFn), 


(9) 


where the functions v Ib; : M 15 —> R are nonnegative convex and positively 
one-homogeneous for i = 0 ,,M. Here, F x n : R 3 — > R 3 is for any 
F £ R 3x3 and any n, a £ R 3 defined as (F x n)a := F(n x a). As in [33] . 
we assume for 0 < i < M, VF 6 R 3x3 , Vn £ S 2 

9 i(RF- n) = gi (F,n ) VF £ SO(3), (10) 

9 i ( F ,n) = gi(F, —n), (11) 

As in 01 ], we assume that there is some c > 0 such that 

MA)>c\A\ . (12) 


for all 0 < i < M and all A £ R 15 . We introduce a subspace Q C y x Z of 
functions with “finite interfacial energy”, using a slightly modified version 
of [331 Def. 3.1]. It is given as follows: 

Definition 2.1 (Interfacial energy). For any pair ( y,z ) £ y x Z let Si = 
d*Q{ OH where flj := supp 2 j and d*Qi is the measure-theoretic boundary of 
Hi with outer (measure-theoretic) normal ni. We denote by Q C y x Z the 
set of all pairs (y,z) £ y x Z such that for every 0 < i < M there exists a 
measure := (■) £ _A4(fl;R 15 ) with 


at := n,U 2 Si , H t := Vs,jJ x n l U 2 Si and c* := (cof V Si y)n\ s . . 
The interfacial energy is then defined as 

oo efee. 

Here | Jj | denotes the total variation of the measure Ji. 


FintiViZ) '■— ( 


(13) 


(14) 
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We recall that the function f t j is called interface quasiconvex if 

J fiysv, n)dU 2 > n 2 (T)f(G , m) (15) 

for every surface deformation gradient G, every unit vector m with Gm = 
0 , every planar two-dimensional region T with normal m, every (curved) 
surface S with normal n and every smooth map y : S —» M 3 with bd = 
bdT (where bd<S = bdT denote the relative boundaries of the two two- 
dimensional surfaces) and such that y = Gx for x £ bdT, see [3H 05] • 
A surface energy is called Null-Lagrangian if (|15l) is satisfied with equality. 
Furthermore, it has been shown in |33j that / is an interface Null-Lagrangian 
if and only if / is a linear function of n, Fxn and cof Fn. This motivates the 
definition of interface polyconvexity ([8]) -([5 |) . in the analogy to the definition 
of the standard notion of poly convexity. The set of configurations Q in 
Definition 12.11 is the natural space where an energy of type | 8 ]) ([9]) can be 
defined. Let us remark that the measures H t and q can be expressed as 


/ v d Hi = / Vy (V x v ) dx, / v ■ dc* = / (cof Vy) : Vv dx (16) 

J J J vt J 

for all v £ Gq°(Q; M 3 ). Indeed, for k £ {1,2, 3} and 0 < * < M, we have 

/ [Vy(V x v)] k = / [djy k ejemdev m }dx = / [nid j y k e j (. m Vm\dx 
Jfli JQi J dtti 

= / [Vu] fc j[n x v\jdx = / [(Vy x n)v] k dx , 

J dQi J dili 


since V x Vy = 0. With the notation (cof Vy)jj = bij, we also have 

/ (cof Vy) : (Vv)dx = / [b k jdjV k ]dx = - / [djb kj v k ]dx + / [njb kj v k ]dx 

JJ Ui j£li J d£li 

= / [(cof \7y)n\ k v k dx = / (cof Vy)n-vdx, 

Jan , Jani 


where we used the Piola identity V • (cof Vy) = 0. 

We also note that by the assumption (fl2l) . we have the bound 


ll-^ z II.M(n;R( M + 1 ) x3 ) — CEi nt (y, z). (17) 

for some constant C < oo. Consequently, the norm || 2 || B y(Q ;R M+i) is con¬ 
trolled in terms of the interfacial energy in our setting. On the other hand, 
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the norm ||-D^||_a 4 (q.r(m+i)x 3 ) satisfies the conditions in Definition 12.11 In¬ 
deed, this follows from the choice gi(F,n) = a\F\ = a\F x n| for a > 0. 
Another example of an interfacial energy which is included in the Definition 
m is given by the choice gi(F,n ) = a\ cof Fn\, see [35] for more details. 
Notice that the first example penalizes surface gradients which are noncon¬ 
stant along interfaces while the latter one increases with the area of the 
interface. 


Body and surface loads: We assume that the body is exposed to possible 
body and surface loads, and that it is elastically supported on a part To of 
its boundary. The part of the energy related to this loading is given by a 
functional L £ C' 1 ([0, T]; VD 1,P (D; R 3 )) in the form 


L(t, y) := [ b(t ) ■ y dx + [ s(t ) • y d U 2 (x) + 

Jn J ri z 


I y - yD(t)\ 2 <m 2 (x). 

(18) 


Here, 6(t, •) : Q —>• M 3 represents the volume density of some given external 
body forces and s(t, ■) : Ti C dfl —> M 3 describes the density of surface forces 
applied on a part Ti of the boundary. The last term in (fTHl) with yo(t, •) £ 
IT 1,P (D;]R 3 ) represents energy of a spring with a spring stiffness constant 
K > 0. Thus our specimen is elastically supported on To in such a way, that 
for K —> oo y is forced to be close to yr> on To in the sense of the L 2 (To;M 3 ) 
norm. A term of this type already appeared in [35] and its static version also 
in [38]. Namely, prescribing a boundary condition from IT 1_1 / p,p ((9D; M 3 ) 
[36] . it is generally not known whether it can be extended to the whole D 
in such a way that the extension lives in y. It is, to our best knowledge, an 
unsolved problem in three dimensions and therefore it is generically assumed 
in nonlinear elasticity that such an extension exists; cf. m, for instance. 
The last term in (1181) overcomes this drawback. Namely, if yn cannot be 
extended from the boundary as an orientation-preserving map the term in 
question will never be zero regardless values of K > 0. 


2.2 Dissipation 

Evolution is typically connected with dissipation of energy. Experimental 
evidence shows that it is a reasonable approximation in a wide range of rates 
of external loads to anticipate a rate-independent dissipation mechanism. In 
order to set up such a process, we need to define a suitable dissipation func¬ 
tion. Since we consider rate-independent processes, this dissipation will be 
positively one-homogeneous. We associate the dissipation to the magnitude 
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of the time derivative of z, i.e., to |£|m+1j where | • \m+i is a norm on M m+1 . 
Therefore, the specific dissipated energy associated to a change of the variant 
distribution from z 1 to z 2 is postulated as in [19] 

D{z 1 ,z 2 ) := jz 1 - z 2 \ m +i- (19) 

Then the total dissipation reads 

V^z 1 ^ 2 ) := f D(z 1 (x), z 2 (x)) dx . 

Jn 

The P-dissipation of a curve z : [0, T] —* BV(El, {0,1}) with [s,t] C [0, T] is 
correspondingly given by (see e.g. [ 22 ]) 

N 

Dissx>(z, [s,t]) ■= sup | ^ V{z(ti- 1 ), z{ti)) : N £ N, s = t 0 <...< tN = t^. 

3 =i 

2.3 Energetic solution 

Suppose, that we look for the time evolution oit y(t ) £ y and t *->• z(t) £ 
Z during a process time interval [0, T] where T > 0 is the time horizon. We 
use the following notion of solution from 122 ], see also [22] HO]: For every 
admissible configuration, we ask the following conditions to be satisfied for 
all t £ [ 0 , T], 

Definition 2.2 (Energetic solution). We say that ( y,z ) £ T x Z is an 
energetic solution to (S,T>) on the time interval [0, T] ift *->• dtE(y(t), z(t)) £ 
L 1 ((0, T)) and if for all t £ [0, T], the stability condition 

£(t,y(t),z(t)) < £(t,y,z) + V(z(t),z) V(y,z)eQ. (20) 

and the condition of energy balance 

£(t,y(t),z(t)) +Diss v(z-, [ 0 ,t]) = So + /q 2§(s,y(s),z(s)) ds ( 21 ) 

where Eq = £(0, y(0), z(0)), are satisfied. 

An important role in the theory of rate-independent solutions is played by 
the so-called stable states defined for each t £ [0;T], We set 

S(t) := {( y,z ) £ y x Z : £{t,y,z) < £(t,y,z) + T>(z,z)V(y,z) £ Q}. 

Note that by (12011 . any energetic solution ( y,z ) is stable for any fixed time. 
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3 Existence of the energetic solution 


A standard way how to prove the existence of an energetic solution is to 
construct time-discrete minimization problems and then to pass to the limit. 
Before we give the existence proof we need some auxiliary results. For given 
N £ N and for 0 < k < N, we define the time increments tk := kT/N. 
Furthermore, we use the abbreviation q := (y, z ) £ Q. Assume that at t = 0 
there is given an initial distribution of phases z° £ Z and y° £ y such that 
q° = 2 °) £ S( 0). For k = 1,..., N, we define a sequence of minimization 

problems 


minimize £(tk,y,z) +V(z,z k 1 ) , (y,z) £ Q . (22) 

We denote a minimizer of (1221) for a given k as ( y k ,z k ) £ Q. The following 
proposition shows that a minimizer always exists if the elastic energy is not 
identically infinite on Q. 

Lemma 3.1. Assume that p > 3, (Ul)- (f7l) . (J9j). (fTTT) - (fT 2 jl hold and let L £ 
C 1 ([0, T]; W 1 ,p (fi; M 3 )). Let := {yo,z 0 ) £ Q satisfy 6(0, y,z) < + 00 . 
Then there exists a solution qjf := (yk,Zk) to (|22l) for each 1 < k < N. 
Moreover, q£ £ S(tk ) for all 1 < k < N. 

Proof. The proof follows the same lines as the proof of [31 .. Thm. 3.3]. We ap¬ 
ply the direct method of the calculus of variations. We denote the elements of 
the minimizing sequence by a lower index in brackets in order to distinguish 
it from the components of z = (z $,..., zm)- Fix k, so that z k ~ l £ Z is given. 
Let {{y(j), - 2 (j))}jeN C Q be a minimizing sequence for £(tk, ■, •) + £>(■, z k_1 ). 
Using the growth conditions ©, (1101) . and in view of the form of L, it follows 
that there is C > 0 such that ||j/(j) llw^fi;® 3 ) + \\ z {j) llBV(f 2 ; M M + 1 ) < C for all 
j £ N. Furthermore, 

sup(||cof Vy (i) || iP / 2 (n . R 3 x 3 ) + II detVy(^ 11 ^/ 3 ^)) < + 00 , 
j 

where p/3 > 1 by our assumption p > 3. Consequently, after taking a 
subsequence, we may assume that y^ —*■ y in W 1 ,p (f 2 ;R 3 ), det Vy^-j —*■ 
detVy in L p ^(£l), cof Vy^) —^ cofVy in M 3x3 ), and z in 

BV(U;R m+1 ). In particular, we have —> z in L 1 (U; {0,1} M+1 ) and 

z £ Z. Moreover, in view of (1161) . ( Jf){j ) converges weakly* in measures 
to as j — > 00 for all 0 < i < M . Standard results for polyconvex 
materials [2j II, [H] show liminfj^oo Eb(y^j\, zr^) > Eb(y,z). Similarly, 
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L(t k ,y(j)) > L(t k ,y) and lim j ^ 00 'D(z^, z k x ) = V(z,z k x ) due 
to the strong convergence of z^ —> z in L 1 (f2;M M+1 ). Finally, 


liminf £) n t {y<j ), Zu\ ) > E int (y,z) 

J —>-00 W/ 

due to [H Thm. 2.38]. Thus, (y, z) £ yxZ. Using weak sequential continuity 
of y i—)■ cofVy and y i-t Vy we see that the limiting measures J* have 
the form of (1131) . This together with a limit passage in the Ciarlet-Necas 
condition (see m Thm. 5]) shows that (y,z) G Q. Namely, y is injective 
almost everywhere in and det Vy > 0 almost everywhere in D. From (1221) . 
one furthermore easily sees that y( v G S(t k ) for all 1 < k < N. □ 

Denoting by B([0, T];T) the set of bounded maps t i-T y(t) G y for all 
t G [0, T], we have the following result showing the existence of an energetic 
solution. 


Theorem 3.2. Let T > 0, p > 3, yc € C 1 ([0, T]; VF 1,P (D; M 3 )), (@}-([7]), 
©, (HD-GZD. Let (y(0),z(0)) G 5(0) and Let there be ( y,z ) G Q such that 
£(0,y,z) < Too. Then there is and energetic solution to the problem (£,T>) 
such that y G B([0,T]-,y), z G BV([0, T]; L 1 ^; M m+1 ) n L°°(0, T; Z). 

Proof. Let q := (y fc , z k ) be the solution of (1221) which exists by Lemma l3Tl 
and let q N : [0, T] —> Q be given by 

N ( t \. = { q k ifte [4,4+i) if k = 0,...,N-l , 

\q% Ht = T. 

Following [22], we get for some C > 0 and for all IV G N the estimates 

lk^llBV( 0 ,T;L 1 (fi;R M + 1 )) — ^ 11^ IIl°°( 0,T;BU(O;R m + 1 )) — C, (24a) 

112/^ llz.°°(0,T’;VlV 1 >3 3 (r2;M 3 )) < C, (24b) 


as well as the following two-sided energy inequality 



q k )d9 < £(tk,q k ) + V(z k , 


z 


k -1 


) 


< [ tk dtm&jde. 

Jtk-l 


£(t k -l,Qk-i) 


(25) 


The second inequality in ([25D follows since is a minimizer of (1221) and by 
comparison of its energy with q := q k _i- The lower estimate is implied by 
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the stability of qj?_i E §(£&_i) when compared with q : = q]f. Having this 
inequality, the a-priori estimates and a generalized Helly’s selection principle 
HB Cor. 2.8] we get that there is indeed an energetic solution obtained as 
a limit for N —> oo. In particular, the fact that det Vy > 0 a.e. in follows 
from the fact that if tj t, (y(j),Z(j)) £ §(£_/) and (y{j),z^) — k (y,z) in 
W 1,p (fi;M 3 ) x then (y,z) G S(t). Indeed, in particular we 

have Z(j S — > z in L 1 (H;M A/+1 ) and hence for all (y,z) E Q, we get 

£(t,y,z ) < liminf £(t«, yi~\, Z(~\) < liminf £(tj, y, z) + liminf z) 

j—>oo j—>oo j—>oo 

= £(t,y,z)+V(z,z) . 

In particular, as £(tj,y, z) is finite for some (y, z) G Q we get £(t, y, z ) < +oo 
and thus det Vy > 0 a.e. in H in view of Q. □ 

Remark 3.3. Adding a term of the form F i—>• | cof .F| p /(det F) p ~ l , which 
is polyconvex, to the bulk stored energy density we can even show injectivity 
of deformations everywhere in for all time instants. See e.g. |?j Rem. 1.2] 
and also m where such term already appeared. 
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